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Quantum states of the diagonal Bianchi type IX model with negative cosmological constant A 
are obtained by transforming the Chern-Simons solution in Ashtekar's variables to the metric rep- 
resentation. We apply our method developed earlier for A > and obtain five linearly independent 
solutions by using the complete set of topologically inequivalent integration contours in the required 
' generalized Fourier-transformation. A caustic in minisuperspace separates two Euclidean regimes 

. at small and large values of the scale parameter from a single classically interpretable Lorentzian 

0^ ■ regime in between, corresponding to the fact that classically these model-Universes recollapse. Just 

one particular solution out of the five we find gives a normalizable probability distribution on both 
branches of the caustic. However, in contrast to the case of positive cosmological constant, this 
■ particular solution neither satisfies the semi-classical no-boundary condition, nor does the special 

initial condition it picks out for h —> evolve into a classically interpretable Universe. 



o 



I. INTRODUCTION 



. Quantum gravity with a non- vanishing cosmological constant formulated in Ashtekar's spin-connection variables 
' jl|~f| has interesting physical states given by the exponential of the Chern-Simons functional |]-[7j] and appropriate 
transformations thereof. In order to elucidate the physical meaning of such states it is interesting to consider their 
restrictions to spatially homogeneous cosmological models. In a recent paper ||, henceforth quoted as [I], we considered 
the diagonal Bianchi IX models with positive cosmological constant from this point of view and found that five linearly 
independent physical states in the metric representation could be derived from the Chern-Simons functional. This set 
of solutions was found to be in one to one correspondence with the set of topologically different integration contours 
which exist for the generalized Fourier-transformation from the Ashtekar-representation to the metric representation. 
Due to the positivity of the cosmological constant the quantum states found in [I] describe an expanding (or collapsing) 
I , classically interpretable Lorentzian Universe at large scale parameters. On the other hand, at sufficiently small scale 



parameters the action defined by the exponent of the wavefunction becomes imaginary and can be associated only 
with a quantum mechanically allowed Euclidean Universe. The two "phases" are separated by a caustic surface in 
minisuperspace. It was found that only one of the five linearly independent states defines a normalizable probability 
K^i ' distribution on this caustic, that this state satisfies the no-boundary condition of Hartlc and Hawking |)],|l0| semi- 
classically for h — + (which means on scales large compared to the Planck scale), and that, again for % — ► 0, it picks 
out an initial condition which evolves into a classically interpretable Lorentzian Universe. For details and further 
literature we refer to [I] . 

It is now of interest to consider also what happens for negative cosmological constant, even though it seems very 
unlikely that our Universe has A < (The age-problem resulting from recently measured high values of the Hubble- 
parameter and the measured large age of globular star-clusters and also the observed high density of galaxies with 
large red-shifts seem to call for A > 0). Our motivation is rather to try the method of [I] for a model-Universe which 
recollapses, i.e. for which quantum-mcchanically a classically interpretable Lorentzian evolution phase is bounded 
for small and large values of the scale parameter by Euclidean evolution phases. Therefore, we have to expect the 
appearance of two caustic surfaces in the minsuperspace of these models, one at small and the other at large scale 
parameter. Is there still a wavefunction, or are there even several, which give a normalizable probability distribution 
on these surfaces, and how are these wavefunctions related to the no-boundary state? 

To answer these questions we apply in section II the method of [I] to obtain expressions for again five linearly 
independent physical states and identify the caustic surfaces in minisuperspace. In section III we determine the 
behavior of the absolute square of the wavefunctions on the caustic and identify a single physical state which gives a 
normalizable probability distribution in this way. In section IV we summarize our results. Within the narrow class of 
models we consider here they seem to rule out, with high probability, a classically evolving Universe with A < 0. 
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II. QUANTUM STATES GENERATED BY THE CHERN-SIMONS SOLUTION 



In this section we want to construct solutions of the Wheeler-DeWitt equation for the diagonal Bianchi type IX 
model with a cosmological constant A < 0, 



hd a -<f> , 



+ 3(8^) 2 Ae 6Q }V(a,0±;A) = 0, (2.1) 



where $ := 2tt e 2a Tr e 2/3 and (3 = (fly ) := diag (b+ + V% f3- , 0+ - \/3 0- , -2 



(2.2) 



In this notation d+ and <9_ denote derivat ives with respect to the variables /?+ and /3_, respectively. By writing 
the Wheeler-DeWitt equation in the for m (pTl ) we have assumed a specific factor-ordering, which is su gges ted by a 
supersymmetric extension of the model |12j-|14|. A different factor-ordering is obtained by considering (2.1) with <!> 
replaced by — $. In the present paper, as in [I], we will restrict ourselves to the factor-ordering as in (2.1), while a 
brief comment on the solutions in the second case $ 



-$ is given in appendix A. 



If the expression (2.2) for <f> is inserted into the Wheeler-DeWitt equation (2.1), the following more explicit form is 
obtained 
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(2.3) 



where we have introduced the mean scale factor a := 2 e a . As in the case A > 0, solutions of (2.3) can be obtained by 
a transformation to the Ashtekar representation, where the Chern-Simons functional, restricted to the Bianchi type 
IX case, turns out to be an exact solution. Two of the Fourier integrals which occur in the transformation back to 
the metric representation can be carried out analytically without any loss of generality and afterwards the same one 
dimensional integral representation as in [I] is obtained:]] 



^(k, j3±; X) cx J du exp 



- /(sintt;«,/3±) 



(2.4) 



with ftz; k, 0±) := 2 K 2 e~ 2 P+ z + cosh ( 2 ^ f 3 -) _ z 2 + 2 Ke 2 P+ f _ cos U 2 ^ 

I — z z V 

Here we have introduced the new variable k and parameter A 



ne 



-4/3+ 



k := — A a 2 , 
12 



A — kA 
6n 



(2.5) 



(2.6) 



thereby effectively reducing the number of parameters occuring in (|2.4| ), and we shall also make use of the variables 
Kj defined by 



K, r|" . " K G 



(2.7) 



where the (3j are the entries of the diagonal anisotropy matrix (3. The integration contour C in the integral- 
representation (2.4) can be chosen quite freely, as long as a sufficiently strong fall-off for the integrand and its 
it-derivatives at the borders dC of C is guaranteed. The proportionality factor left open in (|]J) may depend on A and 
will be fixed later. 

While in the case A > the curves of steepest descent of 3tf were of interest, now, due to the different sign of A, 
the curves of steepest ascent lead to suitable integration contours. Moreover, new possibilities for the location of the 
saddle-points occur, which we classify as follows: 



1 Here, in contrast to [I], the total action, including the part which effects the similarity transformation between Ashtekar and 
metric variables, has been denned as the exponent of the integrand. 
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• By choosing |«| sufficiently small at fixed (3±, it is always possible to make all five saddle points of f(z) lie on 
the real axis of the complex z-plane, defining the Euclid /-region of minisuperspace. Note, however, that the 
corresponding points in the M-plane of fig|l] (here u = arcsinz) are real-valued only for \z\ < 1, whereas real 
z-values with \z\ > 1 are mapped into complex conjugate pairs of points on the axes Sftu = ±? and periodic 
repetitions thereof. 

• Except for the case (3± = 0, where all five saddle-points are on the real z-axis, there is the possibility for two of 
the saddle-points to become complex in the z-plane, which defines the Lorentzian regime. 

• For large values of |k| one always enters the Euclid //region, where again all five saddle-points of f(z) become 
real- valued. 

Some typical locations of the saddle-points in these different regimes of minisuperspace and the corresponding curves 
of steepest ascent are presented in fig.|l[ 



Euclid I Lorentzian regime Euclid II 




• Saddle-point * Singularity at u=± | — s— <R f increases 

FIG. 1. Saddle-points and curves of steepest ascent of Jft/ in the complex it-plane for A < 0. The picture given in the 
Lorentzian case only holds for K3 = min-ftc,}. The remaining case can easily be constructed by reflecting this figure on the 
imaginary axis. The dashed curves come from —00 with respect to K/ and are given just for completeness. 



By passing from one of these regions to another, a marginal situation occurs, where two of the saddle-points confluate. 
We will refer to the corresponding hypersurface in minisuperspace as the caustic; it has been calculated and is plotted 
in fig.|^. In contrast to the case A > the caustic obtained here consists of an upper and a lower branch, which are 
connected just by a single point at k — —2, (3± — 0. Furthermore, there are kinks at /?+ > 0, /3_ = and also at the 
other half-rays of the /3±-plane, related to the former by the typical /3±-symmetries of diagonal Bianchi IX. 

Obviously, an exactly isotropic Universe f3± — has to stay purely Euclidean throughout its evolution. On the 
other hand, "large" Universes with Lorentzian geometry must become very anisotropic. Apart from the possibility of 
a negative cosmological constant very close to zero, which would allow for large scale parameters even at |«;|-values of 
order one, it seems impossible for the model under investigation to describe the Universe observed today. 
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FIG. 2. The caustic in minisuperspace for A < 



Nevertheless, let us now construct a basis of solutions to t he W heeler-DeWitt equation (2.3) by choosing topologically 
independent integration contours C in the representation (2.4). 




FIG. 3. Basis set of integration curves 
Using the curves defined in fig.@ we introduce the following solutions 
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ge{-,+,3,NB} , (2. 



with ^ := ^ , 



(2.9) 



which, by definition, are real-valued and normalized in accordance with 
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* e (a = 0) = l , ge{0,-,+,3,NB} 



(2.10) 



The functions Kq and Jo occuring in (2.8) are the usual modified Bessel functions with index 0. It will be of some 
advantage to replace the solutions \& + and ^_ by 



tfi := 



, 0- > 
, 0- < 



* 2 := 



*+ , /9- < 

, /3- > 



(2.11) 



Then saddle-point expansions of the integrals (2.S) in the limit A — * 0, a and 0± fixed, reveal 



lim * = V° WH 



lim = $° 



A^O 



lim = ie{l,2,3}, 

A^O 



(2.12) 



where the upper index "0" denotes the solutions of the A = 0- model given in [I]. 

Without proof we mention that i e {1, 2, 3}, are three asymmetric solutions which generate each other by cyclic 
permutations of the Kj, so consequently the sum of these states, 



1 3 

(2-13) 



besides and ^nb, turns out to be symmetric with respect to arbitrary Kj-permutations. 

Up to different normalization factors, the asymptotic behavior in the limit k — > — oo can immediately be extracted 
from the corresponding expansions in [I] by taking account of the negative sign of A. The only difficulty lies in the 
determination of the saddle-points, which give the dominating contribution to the different solutions in this limit. A 
detailed calculation finally yields 
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at 0± = 0, i.e. while falls of rapidly for a — > oo, the wavefunction ^jvb is strongly divergent in the same limit. 
Moreover, since ^ NB always gets its dominant contribution from the real saddle-point z > 1 (corresponding to the 
points u NB and u* NB in the complex u-plane of fig-|l| via z — smu), just Euclidean geometries are described by this 
state, so we will reject as a physically relevant solution. Note, however, that it is the only state which satisfies 
the no-boundary condition in the limit h — > 0, a — > 0, hence the name of this wavefunction. 

To give the asymptotic behavior in the limit k — > —oo for the states 4^, ie{l, 2,3}, it will be helpful to consider 

:= \ (*j + **) , Ei jk = 1 , (2.16) 



instead. For these solutions the asymptotic expansions 



h. — J i - 2 KjKk 

7T K t 



exp 



A ' K: 



KjHk 



-ijk 



= 1 



(2.17) 



hold, so they fall off very rapidly for a — > oo (remember the negative sign of A !). 

By considering additional asymptotic expansions for large anisotropy it is possible to show that the four states 
i e {0, 1, 2, 3}, are all normalizable on minisuperspace in the distribution sense (see [I] for a discussion of this point 
for A > 0), i.e. so far we are left with a still four dimensional space of physically interesting solutions. 

However, while in the Lorentzian regime receives saddle-point contributions exclusively from the saddle-points 
at complex z and thus describes a Lorentzian Universe in this part of minisuperspace, the states ^i, ie{l,2,3}, get 
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additional Euclidean contributions of similar order of magnitude from real saddle-points and are therefore hard to 
interpret. 

The classical trajectories which are generated by V^o m the semi-classical limit Ti — > in the Lorentzian regime can 
be computed by solving the equations 

da d9/(z ) d/3± dS/(z ) 

dt da ' df d/3± ' 1 ' 



where we have chosen the lapse-function to be N = |Aa 3 . While the complex saddle-point zq occuring in ( 2.18| ) 
is intended to correspond to the point uq of fig|l], the complex conjugate saddle-point Zq = sinup, which describes 
the time-reversed classical evolution, can be considered with the same right. The corresponding second branch of 
the classical evolution of the Universe is actually needed to define the continuation of a classical trajectory which 
has reached the caustic: in appro aching the caustic the saddle-points Zg and Zq confluate and become real-valued, 



so that, in accordance with ( 2.18 ), the time-derivatives of a and (3± vanish. To continue such a trajectory in time, 
the time- reversed version of ( |2 .18 ) has to be considered. Since the Universe is "reflected" in this way whenever it 
meets the caustic, and since in the generic case the classical trajectories have both of their endpoints on the caustic, 
oscillating Universes are described by ^q. 

The numerical results for the classical trajectories which are obtained in the plane /3_ = of the minisuperspace 
are presented in fig.f§. 



|K| 



* 




-1.5 -1 -0.5 0.5 1 1.5 K + 

FIG. 4. Semi-classical trajectories generated by the complex saddle-points in the Lorentzian regime. For simplicity, we have 
restricted the plot to the plane j3- = 0. The arrows indicate the direction of increasing time t in eq. (2.18). 

We should stress one important difference in the pictures which are obtained for the different signs of (3 + : 

While for /3+ > 0, /?_ = all trajectories run to infinite anisotropy (which is, indeed, a peculiarity of the special 
/3±-direction, corresponding to a kink on the caustic, cf. fig.§), in the case [3+ < the trajectories meet the lower 
branch of the caustic again at a finite /3+-value, representing the general situation. This feature gives rise to the 
existence of a special trajectory with coinciding start- and endpoints, hence describing a Universe that never really 
bec omes Lorentzian. The corresponding points in minisuperspace can be calculated analytically, requiring the solution 
of d2.18|) to be tangential to the caustic, with the result 

4 

k = -v / 20J, f3+ + i(3- = i(ln5-51n2)e 22 T 1 , ne{-l,0,l}. (2.19) 
These points will play an important role in the following section. 



III. BEHAVIOR ON THE CAUSTIC 



Since the classical Lorentzian evolution of the Universe described by the wavefunctions Vt^, i e {0, 1, 2, 3}, is bounded 
by the caustic in minisuperspace, the value of \^\\ on the caustic predicted by the different solutions is of particular 
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interest. In fact |$|^ governs the realization of the different possible histories of the Universe and may thus be 
interpreted as the "initial" value distribution for the classical evolution. 

However, at this stage a new problem arises due to the different branches of the caustic. Since the semi-classical 
trajectories airways can be passed through in both directions, it is impossible to distinguish between their start- and 
cndpoints. The distributions of on the upper and lower branch of the caustic may therefore be considered with 
the same right, and we will always discuss them together in the following. 

The numerical results obtained for \^o\c an d |^s|c on the lower caustic are given in figj|, and figj| shows the 
behavior on the upper caustic, which is very similar for the two different solutions. In the following the additional 
indices "u" and "I" denote the upper and lower branch of the caustic, respectively. 




1-1 1-1 

FIG. 5. Initial value distributions generated by and on the lower caustic(A = —3, Ti = 2n). Like the caustic itself, 
the distributions have kinks in some critical /3±- direct ions, which are partially hidden in the figures. 
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FIG. 6. Initial value distribution of on the upper caustic, normalized to unity at j3± — 0. The numerical plots obtained 
for the different wavefunctions and (and thus for $ defined below) on the upper caustic look very similar, so we restrict 
ourselves to a representation of |$o|c,u- The absolute values taken by l^lc.u at (3± — are given by 1.45-10" 11 , 3.56-10" 13 and 
2.43TCP 11 for the wavefunctions ^o, and ^, respectively (here again A = —3, h = 2ir). Since the lower and upper caustic 
coincide at f3± = 0, it is clear that these values hold for the distributions on the lower caustic, too. That is why we suggest 
to consider the two distributions obtained on the different branches of the caustic as analytical continuations of one another 
through the isotropic point. 
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While on the upper caustic both distributions fall off rapidly with increasing [3± and may be shown to be integrable 
over the /3±-plane, there are /3±-directions on the lower caustic, in which \^\ 2 c i approaches a finite value at infinity. 
Consequently, the wavefunctions on the lower branch are not square-integrable, and hence difficult to interpret as 
probability distributions. Nevertheless, as in the case A > 0, one may construct a new wavefunction as a linear 
combination of the two symmetric wavefunctions and By normalizing and to approach unity in the 
critical /3±-directions, the difference of these new functions is square-integrable on the full caustic. To give an explicit 
expression for the quantum state obtained in this way, we introduce the integrals 



+ T 



da; e 



dx e 



oo oo 

:= / dx sin^sinhi) e ' Icosh2t , £$ 2) (a») := / dx cos(4^sinhi) e ^ cosh2t , 



(3.1) 



which, as far as we know, have no simple representation in terms of tabulated functions. Defining now 

Q(A) := 3ne^ —~ ^ ^- , with „ = J- , 

Ko(-n) 2J (1) (8/i)+ J (2) (8 M ) + e-^4 1) ( M ) 2A 



(3.2) 



the new state can be written in the form 



* := 



l-Q 



(3.3) 



where the overall normalization factor has again been chosen to make f e 1 at a = 0. The behavior of \l/ on the 
caustic has been computed for fig]?]. Taking account of the full distribution, three maxima on the lower branc h o f 
the caustic pick out special initial values for the classical evolution of the Universe. The general representation (2.4) 
easily reveals that the wavefunction becomes arbitrarily sharply concentrated about these maxima in the limit A — > 0, 
i.e. in particular in the limit h — > at fixed A. Consequently, in the semi-classical limit there are just three histories 
of the Universe which occur with significant probability. 
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FIG. 7. Initial value distribution generated by $ on the lower caustic (A 
the upper caustic see fig. 6. 



-3, h = 2ir). For the distribution obtained on 
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In the following we shall be interested in the special points of minisuperspace where the maxima of |\&|^ l arise. Using 
the saddle-point method for A — > the asymptotic behavior of the integrals defined in ( |3.l| ) after some calculation 
yields 

a^o 3 l~2 s , 3 
2^ 



2 - o\/r r (i)(-A)^eA , (3.4) 



and with this result the relation 



* - *o (3.5) 



can be shown to hold at least on the lower caustic. Since is a real-valued, non-vanishing wavefunction, the maxima 
of j^l 2 coincide with the maxima of and using (fTa) they may also be calculated from \&o m the semi-classical 



limit. By performing again a saddle-point expansion for A — > 0, now in the integral represe ntati on (2.8) of Wo, the 
maxima of Wo on the caustic can be calculated analytically with the result given exactly by ( 2.1 9| ) . 



Consequently, within the class of solutions considered here, the only quantum state that is square-integrable on 
the full caustic turns out to predict a Universe which never becomes Lorentzian in the classical limit (albeit classical 
Lorentzian solutions of the Bianchi type IX model with negative cosmological constant actually do exist, cf. figj^)- 

IV. CONCLUSION 

In the present paper we constructed exact quantum states for the diagonal Bianchi type IX model with a negative 
cosmological constant. We found that the method presented in [I] for A > is indeed perfectly applicable to the 
model with A < 0. As for A > it gives five linearly independent solutions, which are generated by the Chern- 
Simons state using topologically different integration contours in the generalized Fourier-transformation to the metric 
representation. Imposing the condition that the wavefunction be normalizable on the caustic, just one wavefunction 
remains, that turns out to have some nice additional properties: It is found to be normalizable in minisuperspace 
in the distribution sense and it respects the symmetries of the Bianchi type IX model. However, this state does not 
satisfy the no-boundary condition in the semi-classical limit in contrast to the case A > 0, and it turns out to predict 
a Universe that never becomes Lorentzian, after all. Hence we obtain the result, that, if one allows for a non-zero 
cosmological constant at all, it should be positive, at least as far as the Chern-Simons functional related states of the 
quantized Bianchi IX model are concerned. 
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APPENDIX A: SOLUTIONS IN A DIFFERENT FACTOR-ORDERING 

For completeness, and in order to obtain an important argument for the factor-ordering chosen for the Wheeler- 



DeWitt equation (2.3), we shall make some comments on a further class of solutions, which can again be discussed 
by using the methods of [I]. 

Considering the Wheeler-DeWitt equation in the form ( [2.l| ) one may ask, why we have not chosen a different 
factor-ordering, which is obtained by changing $ — > — $. Th is c hoice, of course, would not have affected the classical 
Hamiltonian, but the quantum correction — — a 2 Tre 2 ^ in ( |2.3| ) would have changed its signj^] Since the coordinate 
transformation a — > i a. A — > —A has exactly the same effect as the above mentioned change of the factor-ordering, 
it is possible to discuss the solutions of the Wheeler-DeWitt equation in the new factor-ordering by considering still 



2 The factor-ordering obtained in this way corresponds to the .4 + -representation introduced by Kodama in ||, in contrast to 
the ^"-representation, which we have considered up to now. 



9 



equation (2.3), but substituting formally a — > ia, A— * —A i n th e solutions. In the following it will be more convenient 
to use the coordinates Kj and A introduced in (2J3) and (2/7), which transform like Kj — > Kj,X — > —A under this 
substitution. 

It should be clear that the solutions of the Whccler-DeWitt equation (2.3) are still of the form (2.4), but while we 
looked at the cases k > 0, A > in [I] and k < 0, A < in the present paper, now the remaining sectors n > 0, A < 
and k < 0, A > are of interest, which, because of the formal substitution A — > —A just mentioned, describe solutions 
for positive and negative cosmological constant in the new factor-ordering, respectively. It is easily che cked that the 
location of the saddle-points, and therefore the caustic, depends only on J(z;k,(3±) defined in eq. ( |2.5| ). This means 
that, irrespective of the sign of A, we deal with the caustic of [I] in the case k > 0, and with the caustic fig.|| in the 
case k < 0. On the other hand, it is just the sign of A which decides whether the integration curves of [I] (for A > 0) 
or of fig.| (for A < 0) give suitable integration contours. 

However, constructing the solutions for the new factor-ordering in this manner and applying the saddle-point 
method to the integral representation fl2.4| ) in the limit of large anisotropy (3±, it finally turns out, that any solution 
to the Wheeler-DeWitt equation in the new factor-ordering diverges for (5± — > oo, at least in some /3±-sectors. In 
other words, in the new factor-ordering there is no solution which is normalizable in minisuperspace, not even in the 
distribution sense. Furthermore, if the behavior of the wavefunctions on the caustic is considered, actually none of 
these solutions is found to be square- integrable with respect to f3± . 

Comparing these result s wi th the nice normalizability properties of the solutions of the Wheeler-DeWitt equation 
in the factor-ordering of (2.1) presented in [I] and the present paper, we believe to have a compelling argument to 
rule out the new factor-ordering. It would be interesting if this argument could be extended even to the general, 
inhomogeneous case of quantum gravity. 
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